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Abstract . Fourteen-year time series of mesospheric and upper stratospheric temperatures from 
the Halogen Occultation Experiment (HALOE) are analyzed and reported. The data have been 
binned according to ten-degree wide latitude zones from 40S to 40N and at 10 altitudes from 43 
to 80 km — a total of 90 separate time series. Multiple linear regression (MLR) analysis 
techniques have been applied to those time series. This study focuses on resolving their 1 1-yr, 
solar cycle (or SC-like) responses and their linear trend terms. Findings for T(z) from HALOE 
are compared directly with published results from ground-based Rayleigh lidar and rocketsonde 
measurements. SC-like responses from HALOE compare well with those from lidar station data 
at low latitudes. The cooling trends from HALOE also agree reasonably well with those from 
the lidar data for the concurrent decade. Cooling trends of the lower mesosphere from HALOE 
are not as large as those from rocketsondes and from lidar station time series of the previous two 
decades, presumably because the changes in the upper stratospheric ozone were near zero during 
the FIALOE time period and did not affect those trends. 
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32 1. Background and Objectives 

33 Earth’s upper atmosphere has been cooling and its pressure surfaces contracting due to the full 

34 radiative effect of the rising greenhouse gas concentrations, particularly of CO 2 [Lastovicka et 

35 al, 2006]. Model predictions are qualitatively consistent with the observed contraction of the 

36 thermosphere based on several decades of satellite drag measurements [Roble and Dickinson, 

37 1989; Emmert et al, 2004]. Remsberg [2008a] also found a long-term cooling on constant 

38 pressure surfaces of the upper mesosphere that was somewhat greater than -1.0 K/decade at 

39 middle latitudes, based on his analyses of 14-yr (1991-2005) temperature time series from the 

40 Halogen Occultation Experiment (HALOE) [Russell et al, 1993] of the Upper Atmosphere 

41 Research Satellite (UARS). That rate of cooling for the mesosphere is also roughly consistent 

42 with what was predicted [ Roble and Dickinson, 1989]. 

43 

44 In the present paper the temperaUire trends from the HALOE data are reanalyzed and reported 

45 for constant altitudes rather than pressure levels, in order to compare them more directly with 

46 some recent published findings from ground-based lidar and microwave measurements for the 

47 same time period. Trends in temperature versus altitude (or T(z)) are affected by temperature 

48 changes in the underlying atmospheric column, while the results for temperature versus pressure 

49 (or T (p)) are based on physical processes at the local pressure level. The concurrent and 

50 periodic, 1 1-yr response in T(z) is resolved, too, and shown to be in-phase with the solar cycle 

51 flux at most altitudes. The trends in T(z) from the HALOE data are also contrasted with 

52 published results based on lidar and rocket measurements of earlier decades, when decreases in 
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53 ozone were occurring in the stratosphere and contributing to a further contraction of the 

54 atmospheric column from the upper stratosphere to the middle mesosphere. 

55 

56 2. Data Analysis Approach 

57 Remsberg [2008a, and references therein] concluded that the atmospheric sampling of the 

58 HALOE experiment was adequate for resolving the seasonal and longer-period variations in its 

59 temperature profiles after they had been bin-averaged within latitude zones. The 95,900 sunrise 

60 (SR) plus sunset (SS) scan profiles have good signal-to-noise for the altitude range of this study, 

61 and there is a calibration against an exo-atmospheric. Sun-look for each scan [Russell et al., 

62 1993]. Furthermore, no significant trends were found for the HALOE instrument and its 2.8-pm 

63 CO 2 channel transmission measurements that were used for the retrieval of T(z) above about 38 

64 km [ Gordley et al, 2009]. Below that altitude the retrieved HALOE T(z) essentially relies on a 

65 tie-on to profiles from the 12Z operational temperature analyses provided to the UARS Project 

66 by the NOAA Climate Prediction Center (CPC). A retrieval tie-on was also made to the Mass 

67 Spectrometer and Incoherent Scatter (MSIS)-90 empirical atmospheric model above about 85 

68 km. The vertical resolution of the retrieved portion of the HALOE T(z) profile is about 3.5 km 

69 or similar to that reported for ground-based lidar profiles. 

70 

71 Time series of bin- averaged profiles were obtained from the SR and then the SS measurements 

72 within 10-degree wide latitude zones from 40S to 40N and for 10 altitude levels from 43 to 80 

73 km. A minimum of 5 profiles was required for each bin-averaged point. The seasonal sampling 

74 from HALOE was not as good for higher latitudes, so those regions were not evaluated for their 
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75 temperature trends. The SR and SS data at a latitude and altitude were combined into a single 

76 time series after adjusting the SR and SS points for the average difference of their separate time 

77 series. The effects of the diurnal temperature tide are accounted for by that adjustment to first 

78 order, as shown in Remsberg [2007; 2008a]. As an example. Figure 1 presents the adjusted T(z) 

79 time series of over 200 points from HALOE for 30N and 75 km. It is realized that the tidal 

80 amplitudes also undergo some seasonal variation, especially at low latitudes. Thus, the present 

81 approach to a tidal adjustment imparts a small bias to the amplitudes of the seasonal terms for the 

82 present analyses, but not to those of the longer period terms. 

83 

84 The oscillating curve in Figure 1 is a multiple linear regression (MLR) model fit to the time 

85 series points. The MLR model is composed of constant (Const), annual (AO), semi-annual 

86 (SAO), quasi-biennial (QBO) and sub-biennial (IA) terms, an 1 1-yr or solar cycle-like (SC) term, 

87 and a linear trend (Lin) term. The straight line includes the Lin term that has been plotted 

88 relative to its value midway between 1991 and 2005. The QBO and IA terms have average 

89 periods of 853 days and 640 days, respectively, as determined by Fourier fits to the time series 

90 residuals after accounting for the seasonal terms. Amplitudes of those interannual terms in T(z) 

91 are very similar to their values for T(p) in Remsberg [2008a]. After accounting for the seasonal 

92 and interannual terms the time series of the residuals was checked for any remaining periodic 

93 structure — an important test for the acceptance of the final MLR model. Proxy terms related to 

94 forcings from the eruption of Mt. Pinatubo or the aperiodic El Nino/Southern Oscillation 

95 (ENSO) index were not added to the MLR models because the residuals did not show any clearly 

96 anomalous features. 
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97 


98 The linear trend term from the MLR model of Figure 1 is -2.2 K/decade and is significant at the 

99 99% confidence interval (Cl). This trend value is in reasonable accord with other observations, 

100 as will be shown in Section 4. The response of temperature to the solar cycle has been evaluated 

101 by most researchers by fitting the time series points with a proxy term for the solar flux. Instead, 

102 the approach used in this study is to merely fit a sinusoid of 1 1-yr period and then to check the 

103 phase of its maximum to see whether it was within a year or so of solar cycle flux maximum as 

104 described in more detail in Remsberg [2007], This alternate approach has been used because 

105 there may also be decadal-scale changes in temperature due to dynamical forcings [e.g., 

106 Hampson etal., 2005; Shibata and Kodera, 2005]. Nevertheless as shown in Section 3, the 11- 

107 yr MLR term is in-phase (or very nearly so) for many of the temperature time series. In the case 

108 of Figure 1 the amplitude of the 1 1-yr term is highly significant, and its maxima occur in late 

109 January 1992 and 2003 or 1 .1 years after the approximate times of the maxima for the solar cycle 

110 uv-flux. Some studies are showing that the 1 1-yr solar cycle response in temperature can be 

111 more apparent for a particular season or phase of the QBO cycle [e.g., Matthes et al., 2004], 

112 However, the HALOE dataset provides many fewer points for time series partitioned in those 

113 ways, leading to large reductions in the significance of all their MLR terms. The 1 1-yr and trend 

114 terms of the present analyses are representative of the zonal-mean, annual-average temperature 

115 distribution from HALOE. 

116 

117 3. The 11-yr and SC-like Responses 
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118 A total of 90 separate T(z) time series like that of Figure 1 are analyzed in the present study. 

119 Their MLR models include both an 1 1-yr term and a linear trend term. In this section the 

120 character and significance of their 1 1-yr or SC-like terms are considered. Their quality and 

121 validity are judged based on the magnitudes and phases of the max minus min responses and for 

122 the pattern and continuity of the terms within the altitude/latitude domain. If they are reasonable, 

123 then it is likely that the concurrent linear trend terms of Section 4 are meaningful, too. 

124 

125 Figure 2 is the contour plot of the max minus min responses for the 1 1-yr term. The dark 

126 shading denotes those regions having Cl greater than 90%; lighter gray shadings have a Cl of 

127 between 70 and 90%. Figure 3 is the contour plot of the phase of the 1 1-yr terms. Darker 

128 shading denotes those regions where the 1 1-yr terms have their maximum values within +/-1 .5 

129 years of January 1991 and 2002 — the approximate midrange time for the maxima of the 

130 traditional solar flux proxies. Figure 3 shows that the 1 1-yr term is very nearly in-phase with 

131 solar maxima over most the domain. A notable exception occurs for the region of weak 

132 responses of the mid mesosphere at northern hemisphere middle latitudes (Figure 2), just below 

133 the altitude region having a larger response. In fact, the largest 1 1-yr responses occur in the mid 

134 and upper mesosphere at middle latitudes, and they are in-phase. However, there is some 

135 asymmetry between the southern and northern hemispheres for those responses, quite possibly 

136 due to associated differences from effects of decadal-scale, wave forcings. 

137 
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138 Figure 4 is the contour plot of the solar-cycle (SC-like) max minus min T(z) values that were 

139 obtained by multiplying the 1 1-yr responses of Figure 2 by the cosine of the ratio of the 

140 corresponding phase lead or lag (in years) of Figure 3 to the 1 1-yr period of the SC. In effect, 

141 Figure 4 is the approximate response one would find by performing a regression of the HALOE 

142 T(z) time series against a more traditional solar flux proxy. One can directly compare the T(z) 

143 results of Figure 4 with the analysis results for time series of the ground-based Rayleigh lidar 

144 station measurements. As an example, Table 1 provides the HALOE values at ION and at 20N 

145 along with the lidar results at Gadanki (13.5N) of Sridharan et al. [2008] and at Mauna Loa 

146 (19. 5N) of Li et al. [2008], respectively, and for a similar set of years. At both locations there is 

147 good agreement for the magnitudes of the SC-like responses from HALOE and lidar. On the 

148 other hand, the lidar results at these two latitudes tend to be out-of-phase with the solar cycle in 

149 the upper mesosphere, while the HALOE results are less so. It is noted that the lidar results are 

150 averages of measurements obtained over several nighttime hours, whereas the HALOE results 

151 are strictly from its SR and SS measurements. Perhaps this difference indicates that there is a 

152 decadal or solar cycle effect in the tidal forcings at low latitudes, where the tidal amplitudes are 

153 large. Batista et al. [2008] looked for a solar cycle response in their lidar temperature time series 

154 at 23 S, but they did not find any significant values for the altitudes of 40 to 60 km. Figure 4 

155 shows that the responses near 23 S from HALOE are also no greater than about 0.5 K for those 

156 altitudes. 

157 

158 Table 1 also compares the annually-averaged, SC-like response profile from HALOE at 40N 

159 with that from the lidar at Haute Provence, France (44N), but based on its longer period of 
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160 observations from 1979/2001 [Keckhut et ah, 2005]. Nevertheless, because the SC forcing is 

161 essentially periodic the responses of temperature to it ought to be similar for multiple solar 

162 cycles. Responses from those lidar data were reported separately for summer and winter; 

163 numerical responses were estimated from their Figure 3 and then averaged for its profile entry in 

164 Table 1 . The respective HALOE and lidar response profiles agree well in the upper mesosphere 

165 and near the stratopause, but not in the middle mesosphere or at 43 km where the HALOE- 

166 derived responses are small. It may be that there is more wintertime wave activity affecting the 

167 upper stratospheric and mesospheric column over the lidar site than is the case for the zonally- 

168 averaged responses of HALOE [Hampson etal, 2006]. More specifically, the lidar temperatures 

169 at 44N, 6E are often affected noticeably by sudden stratospheric warming (SSW) events during 

170 wintertime, while the HALOE analyses in the 40±5N bin are somewhat equatorward of that 

171 latitude and less subject to such local, large-amplitude wave forcings. 

172 

173 Figure 3 shows that the 1 1-yr responses of the upper stratosphere from HALOE are in-phase 

174 with the solar flux at the Equator but are more nearly out-of-phase at middle latitudes. This 

175 finding is qualitatively consistent with the weak SC-like responses observed in time series 

176 records of lidar and microwave results from 35 to 45 km and for the middle latitude stations of 

177 the Network for the Detection of Atmospheric Composition Change (NDACC) [ Steinbrecht et 

178 a/., 2009]. The response of temperature to the solar forcing should be nearly the same for T(z) 

179 and T(p), and the contour values and patterns of Figure 4 are very similar to those from HALOE 

180 T(p) in Remsberg [2008a, Figure 14]. The SC-like, HALOE responses are of order 1 K in the 

181 region of the tropical stratopause, and the responses within middle atmosphere. 
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182 chemistry/climate models are in good agreement with that value [e.g., Austin et al., 2008; Rind et 

183 al, 2008; Tsutsui et al, 2009]. 

184 

185 4. Trends in T(z) 

186 Temperature trends in the mesosphere are a direct measure of the radiative effects from the 

187 combined changes in the “greenhouse gas” concentrations, primarily those of C0 2 and 0 3 . 

188 Figure 5 shows the linear trends for ozone on a pressure surface from HALOE, so that one can 

189 know how it was changing for the 1991-2005 period at specific levels of the atmosphere. Those 

190 trends are nearly zero through the upper stratosphere and up to at least 0.2 hPa (~60 km). In the 

191 middle stratosphere the HALOE ozone trends were still decreasing, as discussed in more detail 

192 in Remsberg [2008b]. 

193 

194 The associated linear trends in T(z) are presented in this section and compared with concurrent 

195 findings from station data. Then the HALOE T(z) trends are contrasted with the findings from 

196 lidar and rocket station time series data of the preceding decade or two, when the ozone of the 

197 upper stratosphere was declining at a rate of about -5 to -7%/decade [WMO , 2003], Keeping in 

198 mind that the T(z) trend for a given altitude reflects the effects of the ozone changes within the 

199 underlying atmospheric column, it is likely that there was an enhanced rate of cooling in the 

200 upper stratosphere and lower mesosphere at that earlier time [Akmae\’ et al, 2006]. 

201 
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202 Figure 6 shows the T(z) trends from HALOE in terms of K/decade. In order to be consistent 

203 with the convention for the contours of T(p) in Remsberg [2008a], the solid contours are the 

204 cooling trends, while the dashed contours define the zero and wanning trends. The darker 

205 shading denotes those regions where the confidence intervals (Cl) for the trends are greater than 

206 90%; lighter shading defines regions having Cl of between 70 and 90%. One can see that 

207 significant cooling was occurring for this period across the region of the mid to lower 

208 mesosphere, and that it exhibited good continuity from 40S to 40N. In general, the trends in T(z) 

209 are more negative than those for T(p) of Remsberg [2008a] because the trends for T(z) represent 

210 the effects of the cooling for the underlying atmospheric column. T(z) trends in the upper 

211 stratosphere from the HALOE data vary between 0.0 and -0.5 K/decade and are similar to those 

212 reported for the 35 to 45 km region from the NDAAC station data [Steinbrecht, et al, 2008]. 

213 

214 Table 2 gives HALOE T(z) comparisons with the trend profiles from two low latitude lidar 

215 stations — in Brazil [Batista et al., 2008] and in Gadanki [Sridharan et al., 2008]. Cooling trends 

216 were found over the altitudes of the upper stratosphere and lower mesosphere from all the 

217 datasets, and in most instances the HALOE results agree with the lidar values. For the lidar 

218 comparison in Brazil the difference for the trends at 50 km is outside the combined error 

219 estimates from the two datasets. Still, it is important to remember that the HALOE results are 

220 not specific to the longitude of the lidar measurements of Brazil. 

221 


11 



222 Table 3 compares the trends from HALOE at three altitudes of the mesosphere with the values 

223 for previous decades from rocketsondes at 20N and 30N [Keckhut et al., 1999] and from lidar at 

224 44N [Keckhut et al., 1995]. There is no significant difference for the trends from HALOE and 

225 the comparison datasets in the upper mesosphere (75 km), where CO2 contributes the most to the 

226 effects of the radiative cooling. However, at 55 km the trends from HALOE are significantly 

227 smaller than those of the comparison techniques of the previous decades. Abuaev et al. [2006] 

228 reported globally- averaged trends in T(z) from their radiative model calculations that considered 

229 the changes due to CO2 only (Case 1) and then those due to both O3 and CO2 (Case 2) for the 

230 period 1980-2000. Case 2 cooling trends (-1.5 to -2.0 K/decade) are about a factor of three 

231 greater than for Case 1 (-0.6 K/decade) — differences that qualitatively mimic the ones for 55 km 

232 between the comparison measurements of the earlier decade versus those of HALOE. 

233 

234 5. Summary Findings 

235 The 14-yr HALOE dataset was obtained with a sampling frequency that is adequate for resolving 

236 the seasonal and longer-term variations of T(z) in the upper stratosphere and mesosphere, at least 

237 between the latitudes of about 40S to 40N. Though extending for just over one complete solar 

238 cycle, the HALOE T(z) time series yield 1 1-yr, max minus min temperature responses that are of 

239 order 1 K and in-phase with the solar cycle flux at low latitudes near the stratopause. Similar 

240 responses have been reported from ground-based lidar measurements and from model 

241 simulations. In the upper mesosphere the diagnosed, HALOE T(z) responses are somewhat 

242 larger at middle latitudes than from the models. 

243 
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244 The cooling trends from the HALOE dataset are significant at most latitudes of the middle and 

245 lower mesosphere. They range from -1 K/decade at low latitudes to as much as -2.5 K/decade at 

246 the middle latitudes. Values of order -1 K/decade are reasonably consistent with those reported 

247 from lidar measurements at low latitudes and with those predicted for the radiative effects of the 

248 increasing CCT. On the other hand, the cooling rates diagnosed from HALOE are generally 

249 larger than those from models for the upper mesosphere at the middle latitudes. The HALOE 

250 T(z) trends of the lower mesosphere are smaller than those published from rocketsondes and 

251 lidar measurements of the preceding two decades, presumably because those comparative results 

252 were obtained when the decreasing upper stratosphere ozone was an added factor for the total 

253 radiative cooling response through the lower mesosphere. 

254 
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356 Figure 1. Time series of bin-averaged sunrise (SR, open circles) and sunset (SS, solid circles) 

357 temperatures (in K) from Halogen Occultation Experiment (HALOE) measurements at 30±5N 

358 and 75 km altitude. Terms for the multiple linear regression (MLR) model fit are listed at the 

359 lower left (see text). The oscillating curve is the fit for the complete MLR model, while the 

360 straight line is the fit for just the constant plus linear trend terms. 

361 

362 Figure 2. Contour plot of the max minus min values of T(z) for the 1 1 -yr term. Contour 

363 interval is 0.5 K, and the associated pressure profile (left) is approximate. Terms in regions with 

364 darker shading have confidence intervals (Cl) of greater than 90%, while those of lighter shading 

365 have Cl between 70 and 90%. 

366 

367 Figure 3. Phase of the maximum of the 1 1 -year term (in years) as referenced to January 1991 or 

368 2002. Contour interval is 1 .5 years. Negative values are dashed, and the zero value is dotted. 

369 Terms in the shaded regions have a phase maxima within ±1.5 years of those two dates. 

370 

371 Figure 4. Contour plot of the adjusted, solar cycle (or SC-like) max minus min T(z) values. 

372 Contour interval is 0.5 K, and the zero and negative contours are dashed. 

373 

374 Figure 5. Contour plot of the linear trend terms for HALOE ozone versus pressure (%/decade). 

375 Contour interval is 2%. Negative trends have dashed contours, while the zero and positive trends 

376 are solid. Regions of darker shading have Cl greater than 90%; lighter shading denotes Cl of 

377 between 70 and 90%. 

378 
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379 Figure 6. Contour plot of the diagnosed, linear trend terms for T(z) (in K/decade) from 

380 HALOE. Contour interval is 0.5 K/decade. Negative trends have the solid contours, while the 

381 zero and positive trends are dashed. Regions of darker shading have Cl greater than 90%; lighter 

382 shading denotes Cl of between 70 and 90%. 

383 
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0.4 

-1.5 

0.0 

-0.6 

1.0 

0.0 

1.0 

0.7 

0.5 

0.4 

0.7 

0.4 

0.9 

0.3 

0.6 

0.5 


HALOE, 

40±5N 

OHP, 
(43. 9N), 
79-01 

2.6 

3.0 

2.9 

2.3 

0.6 

1.5 

-0.1 

2.2 

0.4 

2.2 

1.4 

2.2 

1.5 

1.5 

1.0 

0.5 

0.6 

0.0 

0.4 

-2.0 
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389 Table 2 — Comparisons of Concurrent Trends in T(z) in K/decade 


Altitude, km 


HALOE, 

25±5S 

Brazil, 

23S, 

93-06 


HALOE, 

10±5N 

Gadanki, 

13.5N, 

98-08 

60 


-0.7 

-1.4 


-1.2 

-1.9 

53 


-0.6 

— 


-0.8 

-0.3 

50 


-1.0 

-2.3 


-0.7 

-0.2 

47 


-0.5 

— 


-0.2 

-0.5 

43 


-0.4 

— 


-0.4 

0.0 

40 


-0.7 

-1.1 
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391 

392 Table 3 — Comparisons with Past Trends in T(z) in K/decade 


Altitude, 

km 


HALOE 

20±5N 

Rocketsonde 

20N, 

69-91 


HALOE 

30±5N 

Rocketsonde 
3 ON, 

69-91 


HALOE 

45±5N 

Lidar 

44N, 

79-93 

75 


-0.2 

-0.7 


-2.2 

-2.9 


-3.0 

-2.0 

65 


-1.5 

-2.8 


-1.1 

-2.1 


-3.5 

-3.0 

55 


-0.9 

-3.4 


-0.6 

-2.7 


-0.4 

-1.0 
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395 Figure 1 . Time series of bin-averaged sunrise (SR, open circles) and sunset (SS, solid circles) 

396 temperatures (in K) from Halogen Occultation Experiment (HALOE) measurements at 30±5N 

397 and 75 km altitude. Terms for the multiple linear regression (MLR) model fit are listed at the 

398 lower left (see text). The oscillating curve is the fit for the complete MLR model, while the 

399 straight line is the fit for just the constant plus linear trend terms. 
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Max Minus Min T(z) Values for 1 1-yr Term 
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401 

Figure 2. Contour plot of the max minus min values of T(z) for the 1 1 -yr term. Contour 
interval is 0.5 K, and the associated pressure profile (left) is approximate. Terms in 
regions with darker shading have confidence intervals (Cl) of greater than 90%, while 
those of lighter shading have Cl between 70 and 90%. 
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402 

Figure 3. Phase of the maximum of the 1 1-year term (in years) as referenced to January 
1991 or 2002. Contour interval is 1.5 years. Negative values are dashed, and the zero 
value is dotted. Terms in the shaded regions have a phase maxima within ±1.5 years of 
those two dates. 
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Figure 4. Contour plot of the adjusted, solar cycle (or SC-like) max minus min T(z) values. 
Contour interval is 0.5 K, and the zero and negative contours are dashed. 
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Linear Trends for HALOE Ozone (%/decade) 



LATITUDE 


Figure 5. Contour plot of the linear trend terms for HALOE ozone versus pressure 
(%/decade). Contour interval is 2%. Negative trends have dashed contours, while the 
zero and positive trends are solid. Regions of darker shading have Cl greater than 90%; 
lighter shading denotes Cl of between 70 and 90%. 


28 



0.001 


Trend for T(z) in K/decade 



405 

Figure 6. Contour plot of the diagnosed, linear trend terms for T(z) (in K/decade) from 
HALOE. Contour interval is 0.5 K/decade. Negative trends have the solid contours, 
while the zero and positive trends are dashed. Regions of darker shading have Cl greater 
than 90%; lighter shading denotes Cl of between 70 and 90%. 
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